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ABSTRACT: This study in centered on models accounting for stochastic defor- 
mations of sample paths of random walks, embedded either in 1? or in 1? . These 
models are immersed in multi-type particle systems with exclusion. Starting from 
examples, we give necessary and sufficient conditions for the underlying Markov 
processes to he reversible, in which case their invariant measure has a Gibbs form. 
Letting the size of the sample path increase, we hnd the convenient scalings bring- 
ing to light phase transition phenomena. Stable and metastable configurations 
are bound to time-periods of limiting deterministic trajectories which are solution 
of nonlinear differential systems: in the example of the ABC model, a system of 
Lotka-Volterra class is obtained, and the periods involve elliptic, hyper-elliptic or 
more general functions. Lastly, we discuss brieEy the contour of a general approach 
allowing to tackle the transient regime via differential equations of Burgers' type. 

1 Introduction 

We are interested in models describing evolution of sample paths of random walks, 
when they are submitted to random local deformations involving possibly several 
links. Roughly speaking, given a finite sample path, say of size N, forming a not 
necessarily closed curve, the problem will be to characterize the evolution of an 
associated family {Yi,i = 1, . . . , N} of Markov processes in the thermodynamic 
limit as N — > oo. This requires to guess and to find the interesting scalings. 

In a previous study [7J, we considered random walks on a square lattice, deforma- 
tions involved pairs of links and occurred at the epoches of Poisson jump processes 
in continuous time (see section [5] for a more exact definition). The analysis was 
carried out by means of an explicit mapping, which led to view the system as a 
coupling of two exclusion processes. 

Starting from a number of observations, we intend to hint in this paper that the 
model in [7J can indeed be cast into a broader class, the ultimate goal being to 
propose methods of wide applicability concerning the following questions: 

• conditions ensuring Gibbs states and explicit forms of the corresponding 
invariant measures; 

• steady-state equations in the thermodynamic limit as N — > oo, and their so- 
lutions in the case of Gibbs states, but also in situations involving permanent 
currents; 

• hydrodynamic and transient equations, when N is sufficiently large, yielding 
thus a complete picture of the evolution. 

Generalizations of the model in Z 2 can follow two natural trends. First, in modi- 
fying the construction of the random walk. Indeed, in the square, lattice we dealt 
with a 4- letter alphabet. Considering instead a finite alphabet of I letters is then 
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tantamount to constructing random walks with oriented links, whose angular af- 
fixes are multiples of , k = 0, . . . , I — 1. The case I = 2 corresponds to the simple 
exclusion process in Z, and I = 3 yields the so-called ABC model. 

Another possible extension is to relax the constraint that the walk lives in I? and 
to define a stochastic deformation process in higher dimension. In the sequel, we 
shall restrict ourselves to some paradigms in Z 2 and Z 3 . 

In section [2 we define a class of two-dimensional models, together with related 
patterns in Z 3 , in terms of exclusion particle systems. Section |31 is devoted to 
stochastic reversibility of the Markov processes of interest and to the Gibbs form 
of their invariant measure. In section the non-symmetric classical ABC model 
is solved (fundamental scaling, phase transitions, classification of stable configura- 
tions) through the analysis of a Lotka-Volterra differential system. The concluding 
section [5] gives a brief overview of ongoing research about large scale dynamics, 
noncquilibrium and transient regimes. 

2 Model descriptions via exclusion particle sys- 
tems 

Our main objective in this section is to show how the evolution of the sample 
paths of random walks can be fruitfully described by means of particle exclusion 
processes. 

Beforehand, to avoid repetition and clumsy notation, let us emphasize that we 
shall only deal with jump Markov processes in continuous time. So implicitly the 
word transition rate will always refer to some underlying generator. Also, N will 
always stand for the size of the sample path, or equivalently the number of its 
links. 

2.1 Preliminaries 

In Z, the simple exclusion model coincides with the well known KPZ system (see 
e -g- El); which represents a fluctuating and eventually growing interface. This 
system is coded by a sequence of binary variables {tj}, j = 1, . . . , N, depending on 
whether a particle is present or not, with asymmetric jump rates. This system has 
been extensively studied. In particular, the invariant measure has been obtained 
in a closed matrix form solution, for fairly arbitrary parameters and boundary 
conditions 0. Large scale dynamics has also been analyzed ^1], showing Burgers' 
equations [Q. 

2.2 2-dimensional models 

1 ) The triangular lattice and the ABC model. Here the evolution of the random 
walk is restricted to the triangular lattice. Each link (or step) of the walk is either 
1, e 2l7r / 3 or e 4i7r / 3 , and quite naturally will be said to be of type A, B and C, 
respectively. This corresponds to the so-called ABC model, since there is a coding 
by a 3-letter alphabet. The set of transitions (or reactions) is given by 

AB h BA, BC h CB, CA h AC, (2.1) 

p+ q+ r+ 
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where the introduced rates are fixed, but not necessarily equal. Also we impose 
periodic boundary conditions on the sample paths. This model was first introduced 
in 0] in the context of particles with exclusion, and a Gibbs form corresponding 
to reversibility has been found in [2j in some cases. 

2) The square lattice and coupled exclusion model. This model was introduced in 
[7] to analyze stochastic deformations of a walk in the square lattice, and it will 
be referred to from now on as the {r a Tb} model. Assuming links counterclockwise 
oriented, the following transitions can take place. 
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We studied a rotation invariant version of this model when 
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Define the mapping (A, B, C, D) -> (r a , r b ) e {0, l} 2 , such that 



'A - 


.(0,0), 


B - 


-(1,0), 


' C- 


-(1=1), 


D- 


-(o,i). 



Then the dynamics can be formulated in terms of coupled exclusion processes. The 
evolution of the sample path is represented by a Markov process with state space 
the 2N binary random variables {r°} and {r b }, j = 1, . . . ,N, taking value 1 if a 
particle is present and otherwise. The jump rates to the right (+) or to the left 
(— ) are given by 

A£(0 = ^NViA* + r^iA ± + r b r b +ll T + rff|V 1 7 ± , 
(i) = t? f? +1 X± + jf rf +1 A=F + ff rf +l7 ± + rf ff +l7 T . 

Notably, one sees the jump rates of a given sequence are locally conditionally 
defined by the complementary sequence. 

3) An extended stochastic clock model. We propose an extension of the pre- 
ceding model for an arbitrary two-dimensional regular graph. To this end, con- 
sider a random walk composed of oriented links, the affixes of which take values 
cut = exp(2i/c7r/n), k — 1, . . . , n, the n-th roots of unity. There are two different 
situations. 

(a) n — 2p + 1 is odd. Then the walk cannot have a fold of two successive 
links, so that local displacements of edges can only be performed by exchanging 
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two successive links. Let X = {X 1 , . . . ,X n } 7 denote the particle types viewed as 
letters of an alphabet, and let {X kl } be the transition rates. The set of reactions 
is defined by 

x k x l <=♦ X l X k , fee [1, 2p + l], k^l. 

\' k 

These rules provide an extension of the ABC model, which we shall discuss in 
detail in section |3J 

(b) When n = 2p is even, the grammar is altered when two successive links fold, 
so that this elementary transition amounts merely to a rotation of the pleat of 
angle ± — (instead of a tt rotation which would occur when exchanging the two 
links). The situation is thus slightly more complicated and the set of reactions is 
now given by 

X k X l «=► X l X k , k = l,...,n, l^k+p, 

(2.3) 

X k x k+ P _» x k+i x k+ P +i^ k = l,...,n, 

where k +p is taken modulo n. It is worth noting that j k (resp. 8 k ) concern folds 
rotating in the counterclockwise direction (resp. clockwise) and that the number of 
letters of each type is no longer conserved. In other words, odd models give rise to 
pure diffusions with eventual drifts, when even models are truly reaction-diffusion 
models. 



2.3 3-dimensional generalizations of the coupled exclusion 
model on a diamond lattice 

2.3.1 Elementary deformations 

Actually, the diamond lattice formulation of stochastic deformations in Z 3 provides 
several straightforward generalization of the 2-dimensional {r a Tb} model. Indeed, 
between two nodes there are 8 = 2 3 possible links (the jumps of the sample path). 
Let (r a , Tf,, t c ) be the vector of binary components corresponding to a displacement 
in each direction, where a, b and c denote here the three particle families (letters) . 

As in the 2-dimensional model, elementary deformations consist in exchanging 
between neighbouring sites the value of one of the binary components. By con- 
struction this model gives a kind of geometric decoupling between the three types 
of particles. In fact all possible existing stochastic coupling result solely from the 
conditional transition rates A^(t;>,t c ), A^(r c ,T a ), and A^("r a , Tf,), according to the 
various possible models. In order to obtain some non-trivial dynamical effect, we 
couple the three systems of particles, keeping in mind the possibility of getting 
back the 2-dimensional {r a Tb} model under certain conditions. This is the subject 
of the next two paragraphs. 
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2.3.2 A linear coupled exclusion model 

Here we propose a coupling which is linear with respect to the complementary 
fields, via the following intensities: 

A±(i) = A±2(7?-7?K, 

A±(i) = A±2(rP-r?)M6, 

A±(t)=A±2(7f-7?)/i c . 

Suppose for a while fi c = and /i a = /i& = /i. Then the sequence {tau c } remains 
disordered, which means the marginal law of t£ is i. As for the subsystem {r a , 17,}, 
up to random contributions (rf — the transition rates A^ and A^ correspond 
to a particular definition of the {r a Tb} model. Hence, in the limit, we obtain a kind 
of disordered {r a Tb} model. 

2.3.3 A non-linear coupled exclusion model 

It is also expected to recover the r a Tb model when one of the components is com- 
pletely frozen in an ordered phase, taking for example r 2 c = for i — 1, . . . , N. 
This situation is fulfilled by choosing the following non-linear couplings 

A±(i) = A±(7? 

Xf(i)=X±(7f 

A±(i)=A±(r? 

2.4 Boundary conditions 

Although they live in Z 2 or Z 3 , let us emphasize that all the objects considered 
throughout this study are curves, hence one-dimensional dynamical systems. 
Hence, for any given sample path of size iV, there are two links referred to as site 
1 and N, at which boundary conditions have to be specified. We shall consider 
only periodic boundary conditions: this means essentially the system is invariant 
under circular permutation of the sites. Consequently, certain geometric quantities 
locally conserved will remain also globally conserved. For example, the distance 
between the two extremities (not necessarily distinct) of a curve remains constant, 
so that closed curves will stay closed for ever. 

3 Reversibility, Gibbs measures and equilibrium 
states 

In the sequel, our goal will be to find exact scalings permitting to derive phase 
transition conditions, for deformation processes of the class defined in section 
as N — > 00. In fact, for the sake of shortness, we shall restrict ourselves to case 
studies where the sample paths have a Gibbs invariant measure, but this is by no 
means a necessary constraint, as commented in section [S] 

Consider a system having a state space of the form W = K s , where if and S 
are finite sets. The following notion, adapted from 1 will be sufficient for our 
purpose. 



-rf)(Tf-ff) Ma , 

-f?M -f?)ti b , 

-rfX^-TfK. 
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Definition 3.1 For any finite set R £ S, let {Vr} be a collection of real numbers. 
A probability measure £ on W is said to be a Gibbs state or a Gibbs measure 
relative to the potential {Vr} if , for all w — (w±, ■ ■ ■ , w\g\) G W, 

R i£R. J 

where Z is a normalizing constant. 

Conditions will now be established, which are either of a geometrical nature or 
bear directly on the transition rates, for the processes of interest to be reversible: 
detailed balance equations hold and they indeed suffice to equilibrate all possible 
cycles in the state space (see e.g. ^0]). In that case a potential does exist and the 
invariant measure can be explicitly expressed as a Gibbs measure. 



CO) = 7j ex P 



3.1 Odd alphabet 

For 2-dimensional walks, when the size n of the alphabet is odd, as already noticed 
in section ?2. 21 there is a dual point of view saying that a n-species particle system 
moves in a one-dimensional lattice: there is exactly one particle per site and the 
transition rates X kl correspond to exchanges of a particle k with a particle / between 
adjacent sites. In a very different context, this model was proposed in [5], from 
which we extract some results pertaining to our topic. 

Up to a slight abuse in the notation, we let X k G {0, 1} denote the binary random 
variable representing the occupation of site i by a letter of type k. The state of the 

system is represented by the array X = {X k , i = 1, . . . , N; k = 1, . . . , n} of size 
N x n. Then the invariant measure of the associated Markov process is given by 

P(X) = iexp[-W(X)], (3.1) 

where 

i<j k,l 

Z is the normalizing constant, and 

ykl 

provided that the following condition holds 

]Ta M 7V fc = 0. (3.3) 

Indeed, a typical balance equation reads 

P[...,X* = l,Xi=l,...] x lk 



P[...,X\ = l,X k +1 = 1,...] A« 



exp(a" £ -a kl ), (3.4) 



and relation l|3.3(l proceeds directly from enforcing the above measure to be in- 
variant by translation. 
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3.2 Even alphabet 



When the cardinal of the alphabet is even, say n = 2p, the situation is rendered 
a bit more involved due possible rotations of consecutive folded links. There is no 
longer conservation law for the number of letters in each class, and one should 
instead introduce the quantities 

A k =N k+p -N k , k = l,...,p-l, 

which represent the differences between populations of links with opposite di- 
rections. Moreover, as a rule, some non-trivial cycles in the state-space are not 
balanced (see figure f5~Tf> . unless transition rates satisfy additional conditions. This 
gives rise to the next theorem. 





Figure 3.1: Elementary cycles: fold (a), 3-link motion (b), square loop (c). 



Theorem 3.2 Assume n = 2p and periodic boundary conditions. Then the system 
is reversible if and only if the following conditions are imposed on the rates and 
on the particles numbers: 

i+p-i k 

« n ^t= i - v/ = 1 '---^ 

k=l 
ykl yk-\-p, I 

(Hi) X! A ' l0 S ^Ife = ' fc = ™- 
2#fc+p 
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The result relies on the next lemma. 

Lemma 3.3 In the case of periodic boundary conditions, if the invariant measure 
has a Gibbs form given by and J^3. then the following relationships must 

hold: 



(iv) 



X kl 

a kl - a lk = log -jj^, k = l,...,n, I ^ k+p; 
k+i,k+ P +i_ a k,k+ P = l 0K JL_ k=i , 



(v) there exists a constant a£l such that 
a ki + a k+ P ,i = a ik + a Lk+p = Qj V fc, Z = 1 

p-1 



Proof. We only present the main steps. Condition (iv) in the lemma comes from 
a balance equation of type l|3.4|l . The case I — k+p corresponds to adjacent folded 
links, and, after setting 

U(k,i) = E (a k+1 ' l +a k+P+1J - a kl -a k+ P- l )X l 31 

l j>i+l 

V(k,i) d = f Y,T,^ k+1+al ' k+P+1 ~ <* lk - a l > k+ v)X], 

l j<i 

equation l|3.4|l has to be replaced by 
P[...,X! = 1,X£* = 1,...] 7 k 



p[...,x* +1 = i,x£? +1 = i,...] s k +i 

= exp [a k+p+1 ' k+1 - a k ' k+p -U(k,i)-V(k, *)] , 
which leads to impose condition (v). 

To take into account the invariance by translation, let a denote a circular permu- 
tation er among the indices, such that o~(i) = 1 + i mod (N). Consider 

k ct (x)= Yl Y, aklx ki) xl °(j)> 

l>k,l>N i<j 

the resulting energy obtained after applying permutation a . Then 

N N 



n a (x) = h(X) + x i Y.( aH ~ a ' fc ) iV ' 

k=l 1=1 



where N l is the number of links I. Since N l+P — N l , I = 1, ... ,p is conserved (but 
N is not), the rule (v) leads to the translational invariance in the form of (vi). 
As to reversibility, one can check that (i) and (ii) are necessary to equilibrate the 
cycles depicted in figure I3~T1 Moreover, the cycle condition imposed by a circular 
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permutation is exactly given by (in). Indeed, this cycle is performed by transport- 
ing one particle through the system from site 1 to site N: during this operation, 
a tagged particle X k will encounter the N l particles corresponding to all other 

species / ^ k and the resulting transition weight is then given by 

which in turn, by using (ii), amounts to condition (Hi) after taking the logarithm. 
These three conditions are in fact sufficient to determine the parameters {a } in 
order to solve (iv), (v), (vi), thus ensuring reversibility. 



4 Steady state of the ABC system in the thermo- 
dynamic limit 

As an illustration, we will present a detailed analysis of the thermodynamic equi- 
librium situation in the case of the ABC model. When there are three particle 
species, the form H of l|3.2[) comes to 

H({X}) = a ab AiBj + a hc BiCj + a ca C l A 3l (4.1) 

i<j 

where the constants a ab , a bc , a ca take the values 

a «b = ^ El a bc = log ^_ a ca = log Z_ 

p q r 

and the constraints (|3.3|) now become 

N B a cal N c a ab ' N A a bc ' ^ ' ' 

4.1 Scaling and Lotka-Volterra equations 

In the example of the ABC model [2], we have at hand an explicit analytic expres- 
sion for the invariant measure. In fact, our claim is that equations at steady state 
in the thermodynamic limit can be derived by using a method proposed in [J] for 
the square lattice model, which a priori does not require any explicit knowledge of 
the invariant measure, which is most of the time untractable. 

We shall only sketch the main lines of argument in the case of the ABC model. 
The number of sites N is fixed and we are interested in the steady state behaviour 
as t — » oo. 

For the sake of shortness, let [A, B, C] = {(A i: Cj), % = 1, . . . , N} denote the 
37V-dimensional bolean vector representing the state of the Markov process. Com- 
bining conditionings and couplings, the approach relies on the construction of a 
stochastic iterative rocking- scheme, to generate the global invariant measure of the 
system. 

The algorithm does construct a convergent sequence of random kernels Q^ n \n > 0. 
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51 Set n — and take the system in some fixed state [^4, B, C)^ = [A, B, C]. 

52 At step n, assume we are given a 3-tuple [A, B,C]^ n \ corresponding to an 
admissible configuration (with respect to the ABC model). To construct 
[A, B, C]( n+1 \ one allows each particle family to evolve as an independent 
one-dimensional exclusion processe until reaching equilibrium, with condi- 
tional transition rates compatible with the rates of the original process intro- 
duced in section l2~2l For instance, with regard to type A particles (omitting 
sometimes the indexing by n for the sake of brevity), we have 

(X+(i)=p+B i +r-C i +TB i C i , 
\ + 1) = p-B, + r+a + r B~i~Ci , 

these expressions being similar to the form of l|2.2(l in the {t T{,} model. T is a 
non-zero quantity taking into account the exclusion at site i in the following 
sense: when BiCi = 1, then necessarily A4 = 1. In the free process, letting 
qf denote the random variable representing the conditional probability of 
finding a particle Ai at site i, we have 

$H , with 

U< W AS" (i + 1)' 4 l-qf 

Using the boolean nature of Bi and Ci, we can write 

(n) 

log-gf = a°»cf> -ct*B<?\ (4.3) 

with analogous quantities uj and for B and C species respectively. 

53 Clearly, the equilibrium states thus obtained are not all admissible, since at 
a site i, several particles of different types may coexist. Therefore, we only 
retain configurations satisfying the constraint A[ n+1 ^ -\-B^ n+1 ' +c\ n+1 ^ = 1, 
for all i, the distribution of which, after some algebra, takes the form 

1 N 

Q([(A,B,C)]^ 1 [(AB) C)]W) = ]J^A^ +V W B W +W W C ( 

i—l 

where is a normalizing constant. Then do n <— n + 1 and go to S2. 

Let 7j = it- + + wj"'. Then, by means of transfer and coupling theorems, 
one can prove the existence of the random vectors 

/,,(«) „(») ,,,(«) \ 
(rf, r?, r?) = lim I ^ , i > 1. 

Fundamental scaling We introduce the so-called fundamental scaling defined 

by 

N \NJ' n \n)' n \nJ' 
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where a, /3, 7 are three positive real constant. 

Limiting equations Set from now on x = i/N, for 1 < i < N. Then one can 
show the weak limits 

p a (x) = lim r a xN , p\x) = lim r b xN , p c (x) = lim r% N 

N—>og N—>oc N— >oo 

exist and satisfy the system of deterministic differential equations 

( d P a in \ 

-q^ = Pa{PPc~lPb), 

= Pb{lPa ~ ape), ( 4 - 4 ) 

— = p c (ap b ~ fip a ), 
{ ox 

with the crucial constraints due to the periodic boundary conditions 

p u (x + 1) = p u (x), Vu = a,6, c. (4.5) 

Proof. Starting from l|4.H[l . one applies the law of large numbers and ergodic 
theorems to justify the approximation of finite sums by Riemann integrals, as 
N -> 00. ■ 
ft is amusing to see that l|4.4f) belongs to the class of generalized Lotka-Volterra 
systems. The original Lotka-Volterra model was the simplest model of predator- 
prey interactions, proposed independently by Lotka (1925) and Volterra (1926), 
see for instance ^2]- Nonetheless, in our case the world is less crual and all particle 
types are treated on an equal footing...! 

Now the form of l|4. 4f> lends itself to a reasonably explicit solution in terms of 
special functions. The first step is to remark the existence of two level surfaces 

Pa+Pb + Pc = 1, (4.6) 

P>IpI = «, (4-7) 

where (|4.6() follows at once from l|4.2|l and k is a constant of motion to be deter- 
mined. Using (j4.6(l to eliminate p c , we rewrite Ij4.7|l as 

P a a pl{\~Pa~PbV = «■ (4-8) 

The change of functions u = p a + Pb and v — (3p a — otpb, yields the first order 
nonlinear differential equation 

^ =(l-u)v(u), (4.9) 
ax 

where v(u) satisfies the equation 

(au + v) a (f3u - vf{\ - u y = K{a + 0) a+l3 . (4.10) 
Formally, u(x) can be expressed as 

r u< -^ du 



u(0) (1-U)v(u) 



(4.11) 
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It appears that the ratios in 14.2f> are rational for all finite N, but, since we have 
let N, Na, Nb, — ► oo, they might become arbitary real numbers in the interval 
[0, 1]. When they are rational, (|4.9|) is a polynomial equation, and then u(x] is in 
some sense just a bit more general that an hyperelliptic function, since, after some 
algebra, we are left with integrals of the form 

y [s 2 — a 2 (l — s)~«]ds, 

where p and q stand for positive integers. The particular case a = (3 = 7 is rather 
simple, since then 

— J = (1- U )[ U 2 (1- U )- K ], 
showing « to be a standard Jacobi elliptic function. 

For any given k, the system (|4.4|l admits of a unique solution. In particular, there 
is always a degenerate solution crumbled to the fixed point 

^ a = — — ' P b = — raZ ' ^ c = — I /j 1 , (4-12) 
a + p + 7 a + p + 7 a + p + 7 

and corresponding to the constant 

K ~ ( a + /? + 7 )«+/»+r ( ^ A6) 

The purpose of the next paragraph is to discriminate between solutions of 
and H4.5JI . in order to relate them to admissible limit-points of as N — > 00 

4.2 Stability, fundamental period and phase transition 

To catch a rough qualitative insight into the solution of (|4.4|l , a standard approach 
relies on a linearization of the right-hand side around the fixed point l)4.12[l . This 
yields a linear differential system, whose matrix 



1 



a + /3 + 7 



—0:7 a/3 ' 
07 -af3 
— /?7 0:7 



has three eigenvalues and A± = ±iy / aPj(a + (3 + 7), and the trajectories are 
located on an ellipsoid. But, since the above eigenvalues are purely imaginary, it 
is well known that no conclusion can be drawn as for the original system, which 
might be of a quite different nature. However, it is pleasant to see that the modulus 
of the nonzero eigenvalues plays in fact a crucial role as shown in the next theorem. 

Theorem 4.1 Let s = a + /3 + 7 and n = -. The limit $ = f lim <& N of the ABC 

3 N^oo 

model is deterministic and there is a second order phase transition phenomenon. 
There exits a critical value 

2tt 

oV PaPbRc 

such that if r] > r] c then there are closed non- degenerate trajectories of sa t~ 
isfying with period T(k p ) = | } p € {1, . . . , [jj-]}- The only admissible stable 

$ corresponds 
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- either to the trajectory associated with k\ if rj > r\ c ; 

- or to the degenerate one consisting of the single point {4-l%\/ if V < Vc- 

The proof involves a forest of technicalities and we only sketch the main lines of 
argument. The first step is to switch to polar coordinates 

dof a 

U a = Pa = r COS 9, 

s 

del ~ P . Q 

u b = Pb = rsmv. 

s 

Rewrite (|4.7|l as 

H(r,6) = logK, (4.14) 

with 

H(r, 9) = ct log [V cos + — ] + /31og[rsin0 + — ] + 7log[— — r(cos 6 + sin 6*)] , 

and let r(8, k) be the single root in r of 14.14(1 . Then 6 satisfies the differential 
equation 

$L = G(9,k), (4.15) 

where 

G{8, k) = - [0(a + 7) cos 9 + a[3 sin 2 9 + a{[3 + 7) sin 2 9] 

+ r(9, k) cos 9 sin 6 [{(3 + 2a + 2 7 ) cos 9 + {a + 2/3 + 2 7 ) sin 9] . 
Letting T(n) be the period of the orbit, we have 



T( 



27r d9 



G(0,k) 



The second important step relies on the monotonic behaviour of T(k) with respect 
to the parameter k, yielding the inequality 

T(k) > T(k). 

Observing that 

r(6,K) = 0, V6>e[0,27r], 
we can write, by l(4.15|l . T(k) as a contour integral on the unit circle, namely 

T(k) - -4«s £ z 2^p_ a ) _ 2ia(3] + 2z[2a/3 + 7 (a + /?)] + 7 (/9- a) + 2ia/3' 
or, after a simple calculus, 



T(k) = 2tt. ' 



a/?7 

which leads precisely to the critical value n c announced in the theorem. 
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5 Perspectives 



This paper is the continuation of |7], but is certainly an intermediate step. For 
the sake of shortness, we did restrict ourselves to the thermodynamics of the ABC 
model. Actually, our goal is to analyze the dynamics of random curves evolving in 
Z m (no spatial constraints) or in Z™, when they warp under the action of some 
stochastic deformation grammar. 

In this project will be carried out in the framework of large scale dynamics for 
exclusion processes, and it will mainly address the points listed hereafter. 
More on thermodynamic equilibrium The trick to derive limiting differen- 
tial systems amounts essentially to writing conditional flow equations on suitable 
sample paths, even in the presence of particle currents. These equations involve 
functionals of Markov and they enjoy special features encountered in many sys- 
tems. It might also be interesting to note that most of the Lotka-Volterra equations 
can be explained in the light of the famous urns of Ehrcnfest. 
Phase transition There exists a global interpretation by means of a free energy 
functional with two components: the entropy of the system, and the algebraic 
area enclosed by the curve. It turns out that the contention between these two 
quantities yield, after taking limits lim t — > oolini/v — ► oo (in that order), either 
stretched deterministic curves or Brownian objects when the scaling is of central 
limit- type. 

Transient regime Our claim is that time-dependent behaviour can be treated 
along the same ideas, up to technical subtleties, by means of a numerical scheme 
based on the conservation of particle currents. This should yield a system of Burg- 
ers equations, extending those obtained in |Jj for the symmetric {r a Tt,} model, 
which had the form 
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